Collisional properties of weakly bound heteronuclear dimers 
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We consider collisional properties of weakly bound heteronuclear molecules (dimers) formed in 
a two-species mixture of atoms with a large mass difference. We focus on dimers containing light 
fermionic atoms as they manifest collisional stability due to an effective dimer-dimer repulsion 
originating from the exchange of the light atoms. In order to solve the dimer-dimer scattering 
problem we develop a theoretical approach, which provides a physically transparent and quantitative 
description of this four-atom system in terms of three- and two-body observables. We calculate the 
elastic scattering amplitude and the rates of inelastic processes such as the trimer formation and the 
relaxation of dimers into deeply bound molecular states. Irrespective of whether the heavy atoms 
are bosons or fermions, the inelastic rate can be significantly lower than the rate of elastic collisions. 
Moreover, the measurement of the inelastic rate which is a four-body observable, can be an efficient 
and precise tool for determining three-body observables such as the three-body parameter, positions 
of Efimov states, and their lifetimes. 

PACS numbers: 



I. INTRODUCTION 

For many years the production of ultracold molecules 
was a long term goal of various experiments [U, 0. A 
breakthrough in this field has recently been achieved 
by combining the ability to cool ultracold atomic gases 
to very low temperatures and the ability to associate 
pairs of them into weakly bound dimers by using Fesh- 
bach resonances H i, H i, 0, i, i, 113, [III- This ap- 
proach turns out to be efficient for obtaining homonu- 
clear weakly bound molecules out of both fermionic and 
bosonic atoms. In the former case the Pauli exclusion 
principle for the constituent atoms makes the molecu- 
lar gas remarkably stable [H, [H, , which has allowed 
the observation of Bose-Einstein condensation (BEG) of 
molecules and investigation of their supcrfluid properties 
d, 0, B @- Molecules consisting of bosonic atoms are 
sufficiently stable only if they are very dilute or com- 
pletely isolated from each other, for example, by an opti- 
cal lattice. However, exploration of such molecules, and 
strongly interacting bosonic atoms in general, is still on 
the short list of many experiments motivated by the po- 
tential opportunity to study the Efimov few-body physics 
[Tsl [igI \1T\ . The evidence for the existence of Efimov 
states has recently been obtained in experiments with 
cesium atoms in Innsbruck p^ . IT^ . 

The tunability of the interspecies interaction in het- 
eronuclear mixtures has been demonstrated in a number 
of experiments 0, M, HH, US, IH, Si, HE HI , which is an 
important step towards obtaining ultracold heteronuclear 
molecules. This is believed to be a rewarding project 
owing to a permanent electric dipole moment that they 
posses in their ground state. There exists a wide range 
of proposals to use such molecules for studies of novel 
quantum phases, realizing quantum computing schemes, 



testing fundamental physical concepts, etc. (see, for ex- 
ample, Refs. [H, [3,[iH, and references therein). 

Taking into account the success in associating homonu- 
clear molecules from ultracold atomic samples, it seems 
reasonable to generalize this technique to mixtures of dif- 
ferent atoms. So far, this scheme has been used to create 
weakly bound K-Rb molecules in an optical lattice [ll|. 
The association of atoms into a weakly bound molecular 
state can be applied, in principle, in many other experi- 
ments studying ultracold mixtures with various fermionic 
and bosonic isotopes. The main concern of this scheme is 
the collisional stability of the molecules that are associ- 
ated into highly excited rovibrational states. Stabilizing 
such weakly bound dimers without isolating them from 
each other opens up a number of interesting research di- 
rections such as the study of novel quantum phases in 
bosonic molecular gases [27*1 and in the presently com- 
pletely unexplored fermionic molecular gases. It would 
also shed light on many intriguing questions related to 
the few-body physics with atoms of different masses and 
statistics. 

The theoretical description of elastic and inelastic 
dimer-dimer collisions requires the solution of a four- 
atom problem which, in a general case, can hardly be 
accomplished. An efficient theoretical approach based 
on the short-range character of the interatomic forces al- 
lows the calculation of scattering properties of weakly 
bound homonuclear and heteronuclear bosonic molecules 
consisting of fermionic atoms [13, Eli [3 ■ However, this 
approach is limited to the case where the mass ratio of 
constituent atoms is smaller than a critical value of 13.6 
and the subsystem of one light and two heavy fermions 
does not exhibit the Efimov effect. Then the scattering 
amplitude is a universal function of the atom-atom scat- 
tering length a and the atom masses, and the so-called 
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three-body parameter is not needed. For larger mass ra- 
tios or in the case of bosonic constituent atoms one has 
to take into account the existence of Efimov three-body 
bound states, the three-body parameter comes into play, 
and it is still debated whether an additional four-body 
parameter is required to describe four-body observables 

In this paper we discuss the scattering properties of 
weakly bound dimers which consist of atoms with a large 
mass difference. Both the light and heavy atoms are in 
a single quantum state, and the weakly bound molecu- 
lar states are formed in their gaseous mixture at a large 
and positive scattering length for the light-heavy interac- 
tion. Our attention is focused on dimers containing light 
fcrmionic atoms as these dimers should be coUisionally 
stable due to an effective dimer-dimer repulsion originat- 
ing from the exchange of the light atoms. We develop an 
approach for solving the dimer-dimer scattering problem, 
which in addition to the short-range character of the in- 
teratomic interactions, relies on a large difference in the 
masses of atoms forming a dimer. The latter circum- 
stance allows us to use the Born-Oppenheimer adiabatic 
approximation and integrate out the motion of one light 
atom substituting it by an effective potential acting on 
the heavy atoms. This reduces the dimer-dimer scat- 
tering to a three-body problem, which can be solved by 
developing the method presented in Ref. [30] . 

Such a hybrid Born-Oppenheimer (HBO) approach 
provides a transparent description of the four-body sys- 
tem of two light and two heavy atoms. It allows one to 
express four-body observables in terms of the light-heavy 
scattering length a and the three-body parameter for the 
subsystem of one light and two heavy atoms. Remark- 
ably, this relation can be inverted and can be used to 
retrieve information on the three-body system from the 
collisional properties of the dimers. As we show in this 
paper, the lifetime of a gas of dimers can significantly 
exceed the lifetime of an atomic mixture or an atom- 
dimer mixture. This suggests that indirect measurements 
of three-body effects through four-body observables can 
be more efficient and precise than measurements using 
three-body systems directly. 

We present a detailed analysis of elastic and inelastic 
dimer-dimer collisions, based on the HBO approxima- 
tion. We calculate the s-wave (p-wave) amplitudes of 
elastic scattering between ultracold bosonic (fermionic) 
dimers and discuss the relaxation of dimers into deep 
bound states as well as the formation of trimer bound 
states in dimer-dimer collisions. As shown, the inelas- 
tic rates are significantly suppressed by the long-range 
exchange repulsion between the dimers. 

The paper is organized as follows. In Sec. [TT] we in- 
troduce the Born-Oppenheimer approximation for three- 
and four-body problems in which both light atoms are 
adiabatically eliminated and in Sec. Illll we discuss possi- 
ble inelastic decay mechanisms in dimer-dimer collisions. 
We then justify the need for a more sophisticated HBO 
approach to describe the dimer-dimer scattering. Section 



IIVI contains the derivation of the HBO equation for the 
four-body problem and in Sec. |V] we show how it can be 
used for calculating elastic and inelastic dimer-dimer col- 
lisional properties. We present the results for the dimer- 
dimer scattering amplitude as a function of the atomic 
scattering length and the three-body parameter in the 
case of bosonic and fermionic molecules. We also discuss 
implications of the results for experimental observation 
of dimers, Efimov trimers, and related few-body physics. 
In Sec. IVII we discuss the validity of the HBO approach 
and compare the HBO results with exact calculations for 
moderate mass ratios and in Sec. IVII] we conclude. 



II. BORN-OPPENHEIMER APPROXIMATION 
FOR THREE- AND FOUR-BODY SYSTEMS 

Let us consider the interaction between two weakly 
bound dimers consisting of atoms with very different 
masses and discuss the dimer-dimer scattering proper- 
ties. A natural way to benefit from the large mass ratio is 
to use the Born-Oppenheimer (BO) approximation [sij . 
where one assumes that the state of fast light atoms adi- 
abatically adjusts itself to the positions of the slow heavy 
atoms. In Ref. [33| the BO approach has been developed 
for a three-body system with an emphasis on the Efimov 
effect. The four-body case has been briefly discussed in 
Ref. 

Within the BO approach the four-body problem can 
be split into two parts. First, one calculates the wave- 
functions and binding energies of the light atoms in the 
field of the heavy atoms fixed at a distance R from each 
other. The sum of the light-atom binding energies gives 
the potential energy surface for the heavy atoms, and 
the second step of the BO procedure is to consider the 
motion of the heavy atoms in this effective potential. 

We assume that the interaction between the light 
atoms is not resonant and can be neglected. Then each 
of the light atoms independently moves in the field cre- 
ated by the pair of heavy atoms. For the interaction 
between light and heavy atoms we use the zero-range 
Bethe-Peierls approach |33| assuming that the motion of 
light atoms is free everywhere except for a vanishing dis- 
tance between the light and heavy atoms. This approach 
makes sense for i? ^ i?e and 

a>i?e, (1) 

where a is the scattering length for the light-heavy in- 
teraction and Re, Re are the characteristic radii of the 
light-heavy and heavy-heavy interatomic potentials, re- 
spectively. Then there are two bound states of a light 
atom in the field of two heavy ones: the gerade state 
{+) with the wavefunction remaining unchanged under 
permutation of the heavy atoms (R — R), and the 
ungerade state (— ) with the wavefunction changing its 
sign under this operation. The corresponding wavefunc- 
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tions are given by 

'g-K±(i?)|r-R/2| 



± 



|r-R/2| 



± 



-K±(_R)|r+R/2| 

|r + R/2| 



(2) 



where M± are normaUzation coefRcients which depend on 
R. The corresponding binding energies are 



e±(^) = -K±(i?)/2m, 



(3) 



where m is the mass of a light atom, and we put h = 1. 
The parameters k±{R) follow from the equation 



K±{R) T exp [-K±{R)R] /R = 1/a. 



(4) 



Equation ^ is obtained by using the Bethe-Peierls 
boundary conditions for the wavefunctions ipR of Eq. ^ 
at vanishing light-heavy separations: 



|r±R/2| 



|r±R/2|^0. (5) 



For large R the perturbative expansion of Eq. (g]) up 
to terms of the second order in the small parameter 
exp(— i?/a) leads to the binding energies 

e±{R) « -kol T 2|eo|^exp(-i?/a) + (6) 

where eo = —l/2ma^ is the binding energy of a single 
molecule within the BO approximation with fixed heavy 
atoms and 



U..{R} - 4|eo|^ (l - ^) exp(-2i?/a 



(7) 



The ungerade (— ) state energy is always higher than the 
energy of the gerade (+) state. Moreover, for i? < a the 
ungerade state is no longer bound. The gerade state is 
bound for any R. For i? <C a its energy is equal to 



e+(i?) « -0.16/TOi?2 



(8) 



In the case where the light atoms are noninteracting 
bosons or distinguishable fermions, it is energetically fa- 
vorable for them to occupy the gerade state and follow 
it adiabatically when the heavy atoms modify their rel- 
ative positions. Then the light atoms mediate an effec- 
tive attraction between the heavy atoms. The effective 
potential in this case is 2e+(i?) -f- 2|eo|, and the remain- 
ing part of the BO procedure is to solve the Schrodinger 
equation for the heavy atoms moving in this potential. 
However, such molecules are expected to be very short- 
lived because of the collisional relaxation into deep bound 
states. Qualitatively, this is similar to the situation with 
homonuclear dimers consisting of bosonic atoms. 

From this point on we will focus on the case where light 
atoms are identical fermions, and light-heavy molecules 
can be long-lived ■ Then the wavefunction of two light 
atoms in our BO problem is an antisymmetrized product 
of the gerade and ungerade wavefunctions 

V'R(ri,r2) = [V-^CrOVRlrs) - ^+{r2)^:^{r^)]/V2. (9) 



The BO adiabatic approach is valid at distances R > a, 
where the effective interaction potential between the 
molecules is the sum of e+{R) and e_(i?). At suffi- 
ciently large inter-heavy separations satisfying the con- 
dition exp(— i?/a) <^ 1, the effective potential can be 
written as 

C/eff(i?) = e+{R) + + 2|eo| « U.^iR). (10) 

The potential C/ex originates from the exchange of light 
fermions and thus can be treated as an exchange inter- 
action. It is purely repulsive and, according to Eq. ([7]), 
has the asymptotic shape of a Yukawa potential at large 
R. Direct calculations show that Ucx is a very good ap- 
proximation to UeS for i? > 1.5a. 

We now turn to the second stage of the BO approach 
and consider the relative motion of two molecules in 
the center of mass reference frame. The corresponding 
Schrodinger equation reads 



(-V^/M + C/eff(i?) - e)*(R) = 0, 



(11) 



where e is the collision energy and M is the mass of a 
heavy atom. Note, that the repulsive effective potential 
is inversely proportional to the light mass m, whereas 
the kinetic energy operator in Eq. (|lip has a prefac- 
tor 1/M. Therefore, for a large mass ratio M/m, the 
heavy atoms approach each other at distances smaller 
than a with an exponentially small tunneling probability 
P cx exp{—By/ M/m), where B ^ 1. Our analysis shows 
that the elastic part of the scattering amplitude can be 
calculated with a very high accuracy from Eq. pT|) for 
M/m > 20 and is practically insensitive to the way we 
choose the boundary condition for the wavefunction at 
R = a. 

When the heavy atoms are fermions one has bosonic 
light-heavy molecules and the dimer-dimer s-wave scat- 
tering length add is of the order of a In ^J~Mjm. The ef- 
fective range of the potential has the same property. Let 
us demonstrate the calculation of add for two bosonic 
dimers in the limit M/m 3> 1. In this case the domi- 
nant contribution to the scattering comes from distances 
in the vicinity of R ~ add S> a, where the effective po- 
tential can be approximated by Eq. ([7]) with a constant 
preexponential factor 



C/cff(i?) « 2(maadd)"^exp(-2i?/a). 



(12) 



Then the zero energy solution of Eq. (fTTj) that decays 
towards smaller R reads 



R Xymadd 



(13) 



where Kq is the decaying Bessel function. Comparing the 
result of Eq. at large R with the asymptotic behavior 
^'(i?) (X (1 — add/ R) we obtain an equation for add'- 



a, / e^T M a \ 



add = - In 

2 V 2 m add J 



(14) 
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This gives add ~ a In M/m, and the scattering cross 
section is 

add = &Trald- (15) 

From Eq. (|13p we see that the interval of distances near 
R = Odd, where the wavefunction changes, is of the order 
of a. This justifies the use of Eq. (|12p . In fact, the cor- 
rections to Eq. (|14p can be obtained by treating the dif- 
ference between Eqs. ^ and (|12p perturbatively. In this 
way the first order correction to the dimer-dimer scat- 
tering length is — (3/4) a^/a^d, where add is determined 
from Eq. (fH)) . 

Qualitatively, Ucft{R) can be viewed as a hardcore po- 
tential with the radius add, where the edge is smeared 
out on a lengthscale ~ a ^ add- Therefore, the ultracold 
limit for dimer-dimer collisions, required for the validity 
of Eq. (fT5|) . is realized for relative momenta of the dimcrs, 
fc, satisfying the inequality 

kadd < 1. (16) 

It can be useful (see Ref. jlj]) to approximate the po- 
tential Ucff by a pure hardcore with the radius add- This 
approximation works under the condition fca ^ 1, which 
is less strict than Eq. (flB)) . 

In the case where the heavy atoms are bosons, the di- 
atomic molecules formed by these atoms with the light 
fcrmions are composite fermions and they can scatter 
from each other only with odd orbital angular momenta. 
Then, at low momenta k satisfying the condition of the 
ultracold limit the leading channel is the p-wave 

scattering. At interdimer separations exceeding the ra- 
dius add of the interaction potential Uce{R), but still 
smaller than their de Broglie wavelength the radial 
wavefunction of the relative motion of the dimers takes 
the form ^'(i?) (x k'^{R — SPdd/R"^), where the quantity 
(3dd is the p-wave scattering volume. The scattering am- 
plitude then reads 

/i? - k'Pdd- (17) 

For the hardcorelike potential Ucs{R) one finds that the 
scattering volume is given by: 

Pdd « (l/3)aL, (18) 

with add following from Eq. (|14p . Accordingly, the scat- 
tering cross section is 

al[J « i87r/3)k*a'dd- (19) 

Note that equations p4|) and (fT8)l are obtained in the 
extreme limit M ^ m and assuming the absence of in- 
elastic processes. The accuracy of these results for not 
extremely large ratios M/m and their possible modifica- 
tions due to the presence of inelastic scattering channels 
will be discussed in Sec. [V] 

Let us now mention that in our numerical calculations 
presented in Sec. |V] we do not find resonances in the 



dimer-dimer scattering amplitude, which could appear 
in the presence of a weakly bound state of two dimers. 
Here we give a qualitative explanation of the absence of 
these bound states. Suppose there is such a state with 
energy e — > 0. Then, at distances R > a the wavefunc- 
tion of the heavy atoms should exponentially decay on 
the distance scale ~ a^mjM <C a, since C/eff represents 
a barrier with the height ~ \/ma^. This means that 
the heavy atoms in such a bound state should be local- 
ized mostly at distances smaller than a. The gerade light 
atom is also localized at these distances as seen from the 
shape of the function The motion of the ungerade 
light atom relative to the localized trimer can be viewed 
as scattering with odd values of the angular momentum, 
and due to the centrifugal barrier the bound states of this 
atom with the trimer should be localized at distances ~ a 
from the heavy atoms. In this case one would expect the 
BO approximation to work, since the ungerade light atom 
is moving much faster than the heavy atoms. However, 
this leads to a contradiction, because in the BO approach 
discussed above the ungerade state at interheavy separa- 
tions R < a is unbound. We thus conclude that weakly 
bound states of two dimers are absent. 

Although there are no resonances in the dimer-dimer 
collisions, there are branch-cut singularities in the scat- 
tering amplitude. They are related to the presence of 
inelastic processes and are discussed in the next section. 

III. INELASTIC PROCESSES. QUALITATIVE 
ANALYSIS 

We distinguish two types of inelastic processes in 
dimer-dimer collisions. The first one is the relaxation 
of one of the colliding dimers into a deep bound state, 
the other dimer being dissociated. The second process is 
the formation of bound trimers consisting of two heavy 
and one light atom, the other light atom carrying away 
the released binding energy. We will start with the for- 
mation of trimer states, which in most cases can be called 
Efimov trimers. The existence of such trimer states can 
be seen from the BO picture for two heavy atoms and one 
light atom in the gerade state. Within the BO approach 
the three-body problem reduces to the calculation of the 
relative motion of the heavy atoms in the effective poten- 
tial created by the light atom. For the light atom in the 
gerade state, this potential is e+{R) found in the previous 
section. The Schrodinger equation for the wavefunction 
of the relative motion of the heavy atoms Xui^) reads 

7?x.(R) = [-V^M + e+iR)] x.(R) = e.X.(R)- (20) 

The trimer states arc nothing else than the bound states 
of heavy atoms in the effective potential e+(i?). Accord- 
ingly, they correspond to the discrete part of the spec- 
trum El,, where the symbol v denotes a set containing 
angular (/) and radial (n) quantum numbers. For R a 
the potential e+(i?) is proportional to —1/R^ [see Eq. ([5])] 
and, if this effective attraction overcomes the centrifugal 
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barrier, we arrive at the well known phenomenon of the 
fall of a particle to the center in an attractive po- 
tential. Then, for a given orbital angular momentum /, 
the radial part of Xi, can be written as 

XAR) ^ R^^^^sm{si\nR/ro), i? < a, (21) 

where 

si = v/0.16M/m- {I + 1/2)2. (22) 

The three-body parameter vq fixes the phase of the wave- 
function at small distances and, in principle, depends 
on I. The wavefunction (j21[) has infinitely many nodes, 
which means that in the zero-range approximation there 
are infinitely many trimer states. This is one of the prop- 
erties of three-body systems with resonant interactions 
discovered by Efimov We see that the fall to the 

center is possible in many angular momentum channels, 
provided the mass ratio is sufficiently large. However, 
for practical purposes and for simplicity, it is sufficient 
to consider the case where the Efimov effect occurs only 
for the angular momentum channel with the lowest pos- 
sible I for a given symmetry. This implies that when the 
heavy atoms are fermions and one has odd I, in order to 
confine ourselves to / = 1 we should have the mass ratio 
in the range 14 < M/m < 76. For bosonic heavy atoms 
where I is even, we set ^ = and consider M/m < 39 to 
avoid the Efimov effect for Z > 2. In both cases we need 
a single three-body parameter tq. 

The formation of Efimov trimers in ultracold dimer- 
dimer collisions is energetically allowed only if < 
— 2|eo|. This means that the trimers that we are inter- 
ested in are relatively well bound and their size is smaller 
than a. Therefore, the process of the trimer formation is 
exponentially reduced for large mass ratios as the heavy 
atoms have to tunnel under the repulsive barrier UeS (R) . 
Moreover, this process requires all of the four atoms to 
approach each other at distances smaller than a, and its 
rate decreases with the trimer size because it is more dif- 
ficult for two identical light fermions to be in a smaller 
volume. In order to carry out the quantitative analysis 
of the trimer formation one should properly take into ac- 
count the motion of the ungerade light fermion at i? < a. 
A method which allows one to do this will be presented 
in Sec. |TVl 

From Eq. (pij) one sees that the behavior of the three- 
body system does not change if rg is multiplied by 

A; = exp(7r/s;). (23) 

On the other hand, the dimensional analysis shows that 
the quantity ejy/eo depends only on the ratio a/rp. This 
means that except for a straightforward scaling with a, 
properties of the three-body system do not change when 
a is multiplied or divided by A/. This discrete scaling 
symmetry of a three-body system, which shows itself in 
the log-periodic dependence of three-body observables, 
has yet to be observed experimentally. In the case of 
three identical bosons, where the BO approach does not 



work and one solves the three-body problem exactly [15| , 
the observation of the consequences of the discrete scal- 
ing requires to change a by a factor of A « 22.7, which 
is technically very difficult in ongoing experiments with 
cold atoms. In this respect three-body systems with a 
very large mass difference can be more favorable because 
of smaller values of A. For example, in order to see one 
period of the log-periodic dependence in a Cs-Cs-Li three- 
body system a has to be changed only by a factor of 
A« 5. 

At this point it is worth emphasizing that three-body 
effects can be observed in a gas of light-heavy dimers, 
where the interdimer repulsion originating from the ex- 
change of the light fermions strongly reduces the decay 
rate associated with relaxation of the dimers into deep 
bound states. The trimer formation in dimer-dimer col- 
lisions is very sensitive to the positions and sizes of Efi- 
mov states, and the measurement of the formation rate 
can be used for demonstrating the discrete scaling sym- 
metry of a three-body system. Indeed, this rate should 
have the log-periodic dependence on a and is detectable 
by measuring the lifetime of the gas of dimers. 

Besides the Efimov trimers, one can have "universal" 
trimer states of one light and two heavy atoms, well de- 
scribed in the zero-range approximation without intro- 
ducing the three-body parameter [s^, [11]. In particu- 
lar, they exist at the orbital angular momentum I = 1 
and mass ratios below the critical value, where the Efi- 
mov effect is absent and short-range physics drops out of 
consideration. One of such states emerges at M/m w 8 
and crosses the trimer formation threshold [ttr = — 2|eo|) 
at M/m « 12.7. The existence of this state is already 
seen in the BO picture. It appears as a bound state of 
fermionic heavy atoms in the potential e+{R) for / = 1, 
and in Sec.|V]we present estimates for the formation rate 
in dimer-dimer collisions. The other state exists at M/m 
even closer to the critical mass ratio and never becomes 
sufficiently deeply bound to be formed in cold dimer- 
dimer collisions. The universal trimer states also exist 
for Z > 1 and M/m > 13.6 [35j. However, the trimer 
formation in dimer-dimer collisions at such mass ratios 
is dominated by the contribution of Efimov trimers with 
smaller I. Therefore, below we focus on the formation of 
Efimov trimers. 

Let us now discuss the relaxation of the dimers into 
deep bound states in dimer-dimer collisions. The typical 
size of a deep bound state is of the order of the char- 
acteristic radius of the corresponding interatomic poten- 
tial. We first consider the relaxation channel that re- 
quires one light and two heavy atoms to approach each 
other to distances ~ i?e <C a. Unlike the trimer forma- 
tion, this decay mechanism is a pure three-body process. 
The other light atom is just a spectator. A qualitative 
scenario of this process is the following. With the tun- 
neling probability which is exponentially suppressed for 
large M/m, two dimers approach each other at distances 
R a. Then the heavy atoms are accelerated towards 
each other in the potential e+(i?), and the light atom 
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in the gerade state is always closely bound to the heavy 
ones as is seen from the shape of the function 1/)+. The 
most convenient way to take into account the relaxation 
process that occurs when the heavy atoms (and the ger- 
ade light fermion) are at interatomic separations Re, 
is to consider the three-body parameter rg as a complex 
quantity and introduce the so-called elasticity parameter 
77* — — s;Arg(ro) [17J. As follows from the asymptotic ex- 
pression for the wavefunction (|2ip . a negative argument 
of To ensures that the incoming flux of heavy atoms is 
not smaller than the outgoing one: 

$out/$in = exp[4s,Arg(ro)] = exp[-4?7,] < 1. (24) 

This mimics the loss of atoms at small distances due 
to the relaxation into deep bound states. In the anal- 
ysis of Efimov states, the imaginary part of tq leads to 
the appearance of an imaginary part of e^. This means 
that any Efimov state has a finite lifetime r due to the 
relaxation. For small |Arg(ro)| and for trimer states 
that are localized at distances smaller than a, we get 
T^^/|eiy| — 4|Arg(ro)| — A'q^/si. Strictly speaking this 
fact indicates that it is not possible to separate the re- 
laxation process from the trimer formation because the 
trimers that are formed in dimer-dimer collisions will 
eventually decay due to the relaxation. Nevertheless, as 
we will show in Sec. |Vl both the modulus and argument 
of the three-body parameter can be determined by mea- 
suring the lifetime of a gas of dimers, leading to a num- 
ber of quantitative predictions concerning the structure 
of Efimov states in the three-body subsystem of one light 
and two heavy atoms. 

Another relaxation channel is the one in which two 
light fermions approach a heavy one at distances ~ i?e ^ 
a. This decay mechanism is quite different from the 
ones described above. On the one hand, it is not ex- 
ponentially suppressed for large mass ratios as the heavy 
atoms do not have to approach each other. However, 
there is a suppression of the decay because of the Fermi 
statistics for the light atoms, which strongly reduces the 
probability of having them in a small volume. We estab- 
lish the dependence of the relaxation rate constant on a 
and estimate its dependence on i?e and the atom masses 
without performing a detailed analysis of the short-range 
physics. For a fixed and large inter-heavy separation R, 
the relaxation rate r(i?) is proportional to the proba- 
bility W of finding two light atoms within a sphere of 
radius ~ Re around the heavy atom, multiplied by the 
frequency of relaxation events. This frequency is of the 
order of 1/mR^ since there is no other energy scale at 
distances ^ Re ■ The probability W can be obtained con- 
sidering the fimit ri,r2 — R/2 in Eq. To the 
leading order in exp(— _R/a) we get 

(25) 

and W - |V'R(ri,r2)|'i?^ with h,2 = |ri,2 + R/2| ~ Re 
in the expression for the wavefunction ip-R. This gives an 



estimate for the relaxation rate T{R) as a function of R: 

T{R) - {l/mRl){Rl/a^R'^)exp{^2R/a). (26) 

Finally, r(i?) should be averaged over the motion of the 
heavy atoms, which leads to the relaxation rate constant 
for the heavy- light-light process in dimer-dimcr collisions: 

o„, = /r(«,|*(«,|V«^ij-I=. (27) 

The main contribution to the integral in Eq. (j27p comes 
from distances R^ a where we use Eq. (|T3| for the wave- 
function '^{R). From Eq. ([27|) we see that ahu oc a~^, 
i.e. this relaxation mechanism is strongly suppressed for 
large dimers. 

In the case of fermionic molecules, i.e. p-wavc colli- 
sions, the relaxation rate constant (I27|l should be multi- 
plied by an additional factor {ka)"^ In^ -y/ M/m, where k 
is the scattering momentum. 

Summarizing the results of sections [Til and IIIII we find 
that the simple BO approach works well for calculating 
the elastic part of the dimer-dimer scattering amplitude 
as the main contribution to this quantity comes from 
large distances between the heavy atoms, where the adi- 
abatic approximation for the motion of light atoms is 
valid. This is also the case for the relaxation into deep 
bound states which occurs in the three-body subsystem 
consisting of one heavy and two light atoms. The other 
inelastic processes (trimer formation and the relaxation 
occurring in the system of one light and two heavy atoms) 
require two heavy atoms to approach each other to dis- 
tances R < a, where the adiabatic approximation breaks 
down. Quantitative studies of these inelastic processes 
will provide us with the estimate for the lifetime of the 
gas of dimers. These studies are also needed to establish 
the relation between three- and four-body observables, 
which can be useful for the exploration of three-body ef- 
fects and Efimov physics. In the next section we present 
a modified BO approach to tackle these issues. 

IV. HYBRID BORN-OPPENHEIMER METHOD 

As we argued in Sec. m the BO method works well 
for the gerade fermion. Its state is characterized by the 
wavefunction ^^(r) and energy e+(_R), both adiabati- 
cally adjusting themselves to the motion of the heavy 
atoms. In the usual BO approximation one integrates 
out the light atom by introducing a potential energy sur- 
face (PES) for the heavy atoms, which in this particular 
case is given by e+(i?). One can improve the PES by 
incorporating the so-called adiabatic Born-Oppenheimer 
correction. The idea of this procedure is to retain ^J^{r) 
as the wavefunction of the light atom, but to calculate 
the PES as the expectation value of the full Hamiltonian 
of the three-body system: 

H3 = - V^/A/ - /2^3 + Int.terms. (28) 
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Here fi^ = 2mM/{2M + m), and r is the coordinate 
of the Hght atom relative to the center of mass of the 
heavy atoms. The interaction terms describe the short 
range hght-hcavy interactions equivalent to the boundary 
conditions ([5]). The PES calculated by using ip^{r) of 
Eq. ^ reads: 



(^+1^31^+) 



-4(i?)/2^, 



(29) 



where fj, = mM/{M + m). It is easy to see that in our 
case the adiabatic correction results in the substitution 
m —^ fi in all formulas of the previous sections. In the 
following we will keep this substitution in mind and, in 
particular, use the same notations e±{R) and ep for the 
"renormalized" quantities. 

Strictly speaking, the adiabatic correction does not im- 
prove the accuracy of the BO procedure as we are still 
missing non-adiabatic terms which are of the same or- 
der of magnitude. However, including this correction has 
the advantage that the renormalized quantities e±(-R) 
asymptotically tend to the exact binding energy of a 
dimer, eq = —1/2^0^, which is more convenient for fur- 
ther analysis. 

Once the light atom in the gerade state is integrated 
out, the original four-body problem reduces to a three- 
body problem, which is described by the three-body 
Schrodinger equation 



[H -Vl/2n3- E]'^(R,r) = 0. 



(30) 



Here H is given by Eq. ^U^j, £; = -2|eo| + e is the total 
energy of the four-body system in the center of mass 
reference frame, and e is the dimer-dimer collision energy. 
The interaction of the light atom with the heavy ones is 
included in the form of the boundary condition ([S]) for ^I^, 
and the ungerade symmetry for this atom is taken into 
account by the condition 



*(R,r) = -*(R,-r) 



(31) 



Let us now discuss the symmetry aspects of the re- 
duced three-body wavefunction ^'(R, r). In the case 
where the heavy atoms are identical fermions, we have 
^'(R,r) = -*(-R,r). Combined with Eq. this 
leads to the condition \1/(R, r) = ^'(— R, — r). Therefore, 
^(R, r) describes atom-dimer scattering with even an- 
gular momenta, and for ultracold collisions we are facing 
the s-wave atom-dimer scattering problem. Performing 
a similar analysis for the case where the heavy atoms are 
identical bosons we naturally arrive at the p-wave atom- 
dimer scattering problem. Both cases are consistent with 
our understanding of the dimer-dimer scattering once we 
recall that there exists a "hidden" gerade fermion. 

In order to solve Eq. (|30p we follow the method of Ref. 
[soj . Namely, we introduce an auxiliary function /(R) 



and write down the wavefunction ^'(R, r) in the form: 
^iR,r)^J2 I X.(R)x:(R')i^K„(2r,R')/(R'), (32) 



where 



-K„[r-R72| -K„|r+R72| 

^-(^-'^^)- 4.|r-RV2r4.|r + RV2| ^''^ 



and 



-i^2fi3{E - e^), < E. 



(34) 



For < E trimcr formation in the state v is possible. 
In this case is imaginary and the function (j33p de- 
scribes an outgoing wave of the light atom moving away 
from the trimer. The choice of the sign in Eq. ([M)) en- 
sures that there is no incoming flux in the atom-trimer 
channel. 

One can check directly that the wavefunction ((32| sat- 
isfies Eqs. (pO| and (|3ip for any /(R). It is also easy 
to see that for the case of fermionic heavy atoms, where 
^'(R,r) = -*(-R,r), one has /(R) = /(-R). Then, 
only antisymmetric ^^'6 contributing to the sum in 
Eq. ([5^ . Otherwise, if we choose antisymmetric /(R), it 
gives ^'(R, r) = — 4'(— R, r) and describes bosonic heavy 
atoms with symmetric Xv 

The function /(R) can be obtained by using the Bethe- 
Peierls boundary condition ^ for the wavefunction (j32p . 
Namely, one has to take the limit r ^ R/2 in Eq. ([32]) 
and fix the ratio between the regular term and the co- 
efiicient in front of the singular term proportional to 
l/|r— R/2|. In this way one finds an integral equation for 
/. However, this equation is not convenient for numeri- 
cal calculations, and we modify it by using the following 
procedure. 

We note that taking a finite number of terms in the 
sum over v in Eq. (|32p docs not give any singularities. 
Therefore, the |r — R/2|^^-peak in ^'(R, r) is provided by 
the high-energy part of the sum, and the corresponding 
coefficient does not depend on the shape of the poten- 
tial €+{R) in the Hamiltonian H. Let us introduce the 
function 'I'°(R, r) by the formula 

vI/0(R,r)=^ / xl^(R)x°*(R')A>.(2r,R')/(R'), 

(35) 

where the superscript in the right hand side corresponds 
to the eigenfunctions and eigenenergies of Eq. (^0)1 with 
e+(i?) = eo. The difference *(R,r) - 5'"(R,r) is not 
singular for r — > R/2 and its value in this limit is given 
by 
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L7(R) = *(R,R/2)-«'°(R,R/2)- / ^[x.(R)x:(R')^«.(R,R')-x°(R)x°*(R')^.2(R,R')]/(R')- (36) 

The function ^'"(R, r) can be calculated straightforwardly by turning the sum over v in Eq. (|35p into an integral 
over k: 

5^X°(R)xr(R') - (2vr)-^ / dV<^-^'\ (37) 

and substituting e° ^ k"^ /Ad + eo. In the limit of r ^ R/2 we get 

*0(R, r ^ R/2) = sin2 e{\v - R/2|-i ~ ^/M^^^) f {H) / - Lf{R), (38) 
where 6 — arctan -^71 + 2M/m, and the integral operator L is given by 

L/(R) = P / {g{\R - R'l) [/(R) - /(R')] ± G(v/i?2 + i?'2 - 2R • R' cos 26')/(R')} , (39) 



with the symbol P indicating the Principal Value of the 
integral (see [l3l)- The plus (minus) sign in Eq. ([39| cor- 
responds to the case of fermionic (bosonic) heavy atoms. 
The function G is defined as 



G{X) 



sin26lA/(eo - E)K2{^/M{eo - E)X/ aint 



(40) 

where K2{z) is the exponentially decaying Bessel (Mac- 
donald) function. 

Finally, making a summation of Eq. p6|) with Eq. ((SS)) 
and comparing the result with the boundary condition 
([5]), we find the following equation for /(R): 



L-L' 



47r 



/(R) = o. (41) 



The operators L and L' conserve angular momentum 
and, expanding the function /(R) in spherical harmon- 
ics, we deal with a set of uncoupled ID integral equations 
for each of the radial functions fi{R). These equations 
for / = and / = 1 are given in Appendix. 



V. LOW-ENERGY DIMER-DIMER 
SCATTERING 

The hybrid Born-Oppenheimer approximation (HBO) 
developed in Sec. IIVI can be used in various ways for ob- 
taining the rates of elastic and inelastic dimer-dimer col- 
lisions. For example, substituting the solution of Eq. ([1T|) 
into Eq. (j32p one can, in principle, restore the total wave- 
function ^'(R, r) and find the scattering amplitudes from 
the asymptotic shape of this wavefunction at R —^ oo. 
However, for the problem of ultracold dimer-dimer scat- 
tering this is not needed as there is a more elegant pro- 
cedure. 



At large distances {R 3> a) the reduced wavefunction 
^(R, r) takes the form: 



*(R,r)«*(R)V^^(r). 



Comparing the singular parts of Eqs. 
limit of r ^ R/2 we obtain 



/(R)cx*(R); i?>a. 



and 



(42) 
ti the 



(43) 



Thus, at large distances /(R) can serve as the wave- 
function for the dimer-dimer motion. In particular, it 
contains the dimer-dimer scattering phase shift. 

In the HBO approach the dimer-dimer scattering am- 
plitude add is determined from the large distance be- 
havior of the solution of Eq. (pij) for E = 2eo. In the 
case of fermionic heavy atoms, the leading channel at ul- 
tralow energies is the s-wave scattering, and the solution 
of Eq. ([iTjl should be matched with 



fo{R) cx {l/R-l/add) 



(44) 



at i? ^ a In ^fMjm. In the case of the p-wave dimer- 
dimer scattering (bosonic heavy atoms) one should match 
the solution of Eq. (|^T|) at these distances with 



/i(i?)cx (i?-3/3W-R')■ 



(45) 



The numerical procedure of HBO is described below. The 
dotted curve in Fig. [1] represents a^d/a from Eq. p4p and 
(3/3drf)^/Va extracted from Eq. The HBO results 

for these quantities are shown by the solid and dashed 
curves, respectively. One sees a relatively good agree- 
ment even for moderate values of M/m. The HBO result 
for add agrees with our previous calculations which use 
the exact four-body equation for M/m < 13.6 (see [14 1 
and Sec. IVI|1 . We also get a good agreement with the 
Monte Carlo results for M/m < 20 |3g |. 
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20 40 60 80 100 

M/m 

FIG. 1: The dimer-dimer s-wave scattering length add/a, 
(sohd) and the quantity {ifiddY^^ / a for p-wave scattering 
(dashed) calculated by using HBO. The dotted line is the 
result of Eq. (fT4)) . 

In this formalism it is straightforward to account for 
inelastic processes of trimer formation and the relaxation 
of dimers into deep bound states. In the latter case we 
mean the relaxation process that requires one light and 
two heavy atoms to closely approach each other. The 
heavy-light-light relaxation channel is well described by 
the simple BO approach (see Sec. IIII|) and we do not 
consider it here. 

Let us first assume that the rate of relaxation into deep 
molecular states is negligible and omit this process. Then 
the three-body parameter is real, and the trimer forma- 
tion rate is determined by the imaginary part of the s- 
wave scattering length or the p-wave scattering volume. 
The rate constant is given by [37| 

Q, = X / Iin("d<i)' heavy fermions, 

M 1 3fc^Im(/3(j(i), heavy bosons. ^ ' 

Alternatively, if one needs to know the rate of trimer 
formation in the state v, one can substitute the solution 
of Eq. (gll) into Eq. and calculate the flux of light 
atoms at r — s- cx). Summing over v gives the same result 
as Eq. . We find that the contribution of the highest 
"dangerous" trimer state is by far dominant and a is very 
sensitive to its position. 

We now include the relaxation of the dimers into deep 
bound states. As we have mentioned in Sec. lIIIl the light- 
heavy-heavy relaxation process can be taken into account 
by adding an imaginary part to the three-body parame- 
ter. The total inelastic decay rate is then still given by 
Eq. (|46p . However, strictly speaking we can no longer 
distinguish between the trimer formation in a particu- 
lar state and the relaxation since the trimers ultimately 
decay due to the relaxation process. In this sense the 
only decay channel is the relaxation. However, for a suf- 
ficiently long lifetime of a trimer, i.e. if the trimer states 
are narrow resonances, we can still see a pronounced de- 



pendence of the total inelastic decay rate on the position 
of the highest "dangerous" trimer state (see below). 

The calculation of the intrinsic lifetime of a trimer re- 
quires a detailed knowledge of short-range physics and is 
beyond the scope of this paper. Estimates of the imagi- 
nary part of the trimer energy, r~^, from the experimen- 
tal data on Csa trimers show that it is approximately 
by a factor of 4 smaller than the real part (in this case 
77* ~ 0.06). From a general point of view, we do not ex- 
pect that the trimers with a binding energy t^, < — 2|eo| 
are very long-lived. However, one can have relatively nar- 
row resonances, and we perform calculations for various 
values of the elasticity parameter 77*. 

In the case of a complex three-body parameter the 
Hamiltonian H (j20p combined with the boundary con- 
dition (j2ip . no longer represents a Hermitian operator. 
Instead, this operator is equivalent to the Schrodinger 
operator with a complex absorbing potential, and in a 
real basis set its representation is a complex symmet- 
ric matrix. Methods of calculating Green functions and 
other properties of such operators are known in quantum 
chemistry (see [s^ and refs. therein). The formalism is 
based on using the canonical symmetric bilinear form in- 
stead of the Hermitian inner product. In our case, this 
results in the formal removal of the sign of the complex 
conjugation for all xt in Sec. IIVI This does not lead 
to any contradiction with the "Hermitian" case as the 
eigenvectors of Hermitian operators can be made real. 

We have performed a numerical analysis of the elastic 
and inelastic rates for bosonic and fermionic heteronu- 
clear dimers for different values of M/m. A brief sketch 
of the numerical procedure is the following. We diagonal- 
ize the Hamiltonian H of Eq. ([^0]) on a finite but large 
grid in coordinate space and find the sets {e,y, Xi/(R.)} 
and {e^,x°(R.)}- The latter is obtained from Eq. ^ 
with e+(i?) = 69. The three-body boundary condition 
(I2ip for Xv is implemented by introducing an additional 
(complex) potential at very short distances. Then, we 
calculate the kernels of the integral operators L and L' 
and solve Eq. (|4T|l . In order to characterize the elastic 
and inelastic dimcr-dimer scattering properties we use 
Eqs. (m, (gSl), and dlSl). 

In Fig. [5] we present the results for the inelastic coUi- 
sional rate in the case of bosonic molecules with the mass 
ratio M/m = 28.5 having in mind ^^^Yb-^Li dimers. The 
solid line corresponds to the case of a real three-body pa- 
rameter. It is convenient to introduce a related quantity, 
ao, defined as the value of a at which the energy, e^, of a 
trimer state exactly equals E = 2eo- This new "danger- 
ous" trimer state becomes more deeply bound for a > oo 
and the rate constant steeply grows being proportional to 
the density of states in the outgoing atom-trimer channel. 
The corresponding orbital angular momentum is equal to 
1 and the threshold law reads (see also inset in Fig. [J) 

a cx const + {E - e,,)^/^ oc const + (a - aq)^/^. (47) 

The constant term in Eq. (IT71) describes the contribu- 
tion of more deeply bound states, which is typically very 
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small. In fact, as trimer states become more compact, 
both light atoms should approach the heavy atoms and 
each other at smaller distances where the trimer forma- 
tion takes place. Since they are identical fermions, there 
is a strong suppression of the trimer formation to these 
deeply bound states. 
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FIG. 2: The inelastic rate constant for bosonic dimers with 
M/m = 28.5 as a function of the atom-atom scattering length 
a. The solid line corresponds to the case of a real three-body 
parameter. The results plotted in dashed, dotted, dashed- 
dotted, and dash-dot-dot lines are obtained by taking into 
account the light-heavy-heavy relaxation processes. The val- 
ues of the elasticity parameter r;, =0.1, 0.5, 1, and oo, re- 
spectively (see text). The quantity ao is the value of a at 
which the energy of a trimer state equals E = 2eo and a new 
inelastic channel opens. The inset shows the region a « ao in 
greater detail in order to see the threshold behavior (|47|l . 



The plot is periodic in the logarithmic scale, the mul- 
tiplicative factor being equal to Ai ~ 7.3 consistent with 
Eqs. ^ and ^ (with m ^i). The dashed, dotted, 
and dashed-dotted curves are obtained for ry, =0.1, 0.5, 
and 1, respectively. The corresponding values of the ratio 
'I'out/'I'in are 0.67, 0.14, and 0.02. The flat line represents 
the limiting case of 77, =00 or $out — 0. This case is uni- 
versal in the sense that physical observables depend only 
on the masses and the atomic scattering length. 

For a very weak light-heavy-heavy relaxation, the 
dimer-dimer inelastic collision can be viewed as the for- 
mation of a trimer (with the rate constant a) followed 
by its slow decay due to the relaxation. In this case one 
can think of detecting the trimers spectroscopically. We 
note, however, that even for the conditions correspond- 
ing to the dashed curve in Fig. [21 i.e. for rj^ as small as 
0.1, the decay rate of a trimer w 0.25|ei/| is rather 
fast, which is likely to make its direct detection diffi- 
cult. For larger ry* one cannot separate the formation 
of trimers from their intrinsic relaxational decay, and a 
is practically the relaxation rate constant. Remarkably, 
according to our results, it remains quite sensitive to the 
positions of trimer states (in this case resonances) for val- 
ues of 77* up to 0.5 and even larger. This suggests that 



measuring the lifetime of a gas of dimers as a function of 
a gives important information on three-body observables. 
Moreover, for small 77* one can have a stable molecular 
gas in sufficiently broad regions of a, where "dangerous" 
trimer states are far from the trimer formation threshold. 

We have also calculated the inelastic rate a for other 
mass ratios in the range 20 < M/m < 76. Its dependence 
on the scattering length has the same structure as plotted 
in Fig. [2I The maximum of the rate constant is well fit- 
ted by the formula a^ax 5.8{a/M) exp{-0.87 ^/ M/m). 
The position of the flat line (77, =00) is also well fitted by 
= 1.6(a/M) exp(-0.82^M/?7i). The multiplicative 



factor in the log-periodic dependence is given by Eqs. 
and (1221) with I = 1. 

We employed the same method for estimating the for- 
mation of the universal trimer state with orbital angu- 
lar momentum ^ = 1 at mass ratios M/m > 12.7 but 
below the critical value for the onset of the Efimov ef- 
fect. The rate constant increases with M/m and reaches 
a = 0.2{a/M) on approach to the critical mass ratio. 
This corresponds to the imaginary part of the scattering 
length Inwdd ~ 4 x 10~^a, which is by a factor of 300 
smaller than the real part of add obtained in four-body 
calculations [Ijl . Thus, the formation of this state does 
not change the elastic scattering amplitude add of Ref. 

We can now estimate the rate constants for ^^^Yb- 
^Li dimers with the range Re < 5nm. On the basis of 
Eq. (1271) and the results in Fig. 121 we find that for a = 
20nni the upper bound of the inelastic rate constant is 
ttmax ~ 4 X 10^^'^cm'^/s and it is larger than the rate of 
relaxation into deep bound states in the system of one 
heavy and two light atoms, ahii- For larger a the heavy- 
light-light relaxation can be neglected and the lifetime of 
the gas of dimers is determined solely by the rate constant 
a. The elastic rate constant for a thermal gas with a ~ 
20nm and T ^ lOOnK equals 



V2T/M- 4 X 10-"cmVs. 



(48) 



Here we used the calculated s-wave dimer-dimer scat- 
tering length add ~ 1.4a. We see that aei is much 
larger than a and this inequality becomes even more pro- 
nounced for larger a due to the scaling relations aei oc a^ 
and Qfmax oc a. 

In Fig. 121 we present the results for the inelastic col- 
lisional rate in the case of fermionic dimers with mass 
ratios 14.5 and 22.2, having in mind ^'''Rb-^Li and ^^^Cs- 
^Li dimers respectively. We observe essentially the 
same structure of the graphs as in the case of bosonic 
molecules. However, now we have p-wave scattering and 
the rate constant acquires an additional factor cx (a/c)^, 
where k is the relative momentum of the dimers. The 
multiplicative factor in the log-periodic dependence is 
now given by Eqs. (|23p and (|22p with ^ = 0, since the 
light-heavy-heavy three-body atomic complex exhibits 
the Efimov effect in the configuration with zero orbital 
angular momentum. The dimer-dimer scattering vol- 
umes are equal to Pdd ~ 0.73a'^ for M/m = 14.5 and 
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fidd ~ 0.98a3 for M/m = 22.2. 




FIG. 3: The inelastic rate constant for fermionic dimers with 
the mass ratios M/m = 14.5 and 22.2. Notations are the 
same as in Fig. [2l 

For fermionic molecules with other mass ratios the in- 
elastic rates behave in a similar way as plotted in Fig. [31 
Considering M/m < 39 in order to avoid the Efimov ef- 
fect in the / = 2 channel, the quantities amax and aoo for 
a given relative momentum k can be well fitted by the 
formulas amax = lA.i{a^k'^ /M) exp(-0. 88 v^M/m) and 
aoo = 2.8(a3fc^/M) exp(-0.77\/Af/m). The rate con- 
stant of elastic collisions is given by 

aei « 487r|/3drfpA:VM. (49) 

In the case of Cs-Li dimers we find that for a ^ 
3nm the upper bound Omax greatly exceeds aui- More- 
over, for a thermal molecular gas with a = 20nm at 
T ~ lOOnK, the quantity amax is comparable with 
and both are approximately equal to 4 x 10~"'^'*cm^/s. We 
suggest that in order to optimize the coUisional rates and 
lifetime of a gas of dimers and obtain a^i 3> a one mod- 
ifies the scattering length or temperature, since dimer- 
dimer inelastic and elastic rates strongly depend on these 
parameters. The scaling relations for the rate constants 
read: aei oc a^k^ and amax oc a^k^. 

It is useful to mention that for fermionic dimers the 
relaxation rates in the case of atom-dimer collisions are 
higher than in the case of dimer-dimer collisions. The re- 
laxation rate constant for collisions of Cs-Li dimers with 
Cs atoms is a ^ a/m « 2 x 10~^''cm^/s and it exceeds 
by 3 orders of magnitude the relaxation rate for collisions 



of these dimers with Li atoms. Although these estimates 
are valid within an order of magnitude, they indicate that 
a pure gas of dimers should have a longer lifetime. At 
least, one should avoid mixing Cs-Li dimers with free Cs 
atoms. 



VI. ON THE VALIDITY OF HBO 

Our derivation of the HBO method in Sec. llVl is based 
on two assumptions. First, we rely on the short-range 
character of the interatomic forces. Namely, all our re- 
sults are valid under the condition ([1]). It assumes that 
the range, i?e, of the interatomic van der Waals potential 
is much smaller than the typical de Broglie wavelength 
of atoms, which is of the order of a. In this case we 
can treat the interatomic potential using the zero-range 
Bethe-Peierls approach |33| . 

The second important assumption is the use of the 
Born-Oppenheimcr approximation for integrating out the 
light fermion in the gerade state. We rely on the pres- 
ence of the small parameter m/M. In the beginning 
of Sec. IIVI we explained how the simplest (clamped 
nuclei) approximation can be improved by adding the 
Born-Oppenheimer diagonal correction. The resulting 
so-called adiabatic approximation provides one with the 
"best possible" potential energy surface for the heavy 
atoms. In principle, one can make further improvements 
and systematically expand the wavefunctions and ener- 
gies in powers of m/M. However, it is believed that the 
amount of work that one has to do in order to get the 
next order term in this expansion is comparable to the 
work required for the exact solution of the problem. 

In such a situation it looks optimal to check the validity 
of the HBO results by comparing them with the ones ob- 
tained in exact four-body calculations |J^]. In Fig. [1] we 
compare add for bosonic molecules in the 3D case, calcu- 
lated exactly (solid curve) and by using the HBO method 
(dashed curve). Already for M/m ^ 10 the results agree 
within the limits of accuracy. The exact method requires 
a huge configuration space for M/m « 10 and its ac- 
curacy drops significantly at larger mass ratios. We es- 
timate the relative precision of this "exact" four-body 
calculation at the highest mass ratio in Fig. 0] to be 2 or 
3%, whereas one should have a much higher accuracy in 
order to resolve the trimcr formation for M/m > 12.7. In 
fact, we believe that the HBO method is more accurate 
for M/m > 10. To justify this belief we have performed 
the same analysis in the ID case, where the exact method 
can be pushed to higher mass ratios. The HBO results 
agree very well with the results of the exact methods 
(see inset in Fig. S]) . In the ID case for the considered 
mass ratios there are no trimer states and we neglect any 
inelastic decay. 

In certain aspects the accuracy of the HBO method 
can be limited by the accuracy of the BO approximation 
for the three-body problem. For example, the BO result 
for Ai in the case of Cs-Cs-Li three-body system is 5.3, 
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FIG. 4: Dimer-dimer s-wave scattering length obtained by 
the exact four-body method (sohd hne) and by HBO (dashed 
hne). The inset shows the results of a similar calculation in 
ID. 

whereas the exact value of this multiplicative factor is 
close to 4.9. Accordingly, the results presented in Fig. [3] 
(and also in Fig. ^ can slightly deviate from the actual 
numbers. However, we believe that the HBO approach 
reproduces reasonably well the shape of the inelastic rate 
dependence on a even for these not very large mass ratios 
and, more importantly, it captures the qualitative physics 
behind the processes of trimer formation and relaxation 
into deep bound states. 

Finally, the wavefunction of a four-body system, in 
which one of the three-body subsystems exhibits Efi- 
mov physics, has many nodes (infinite in the zero-range 
approximation) in both real and momentum space. At 
present, as far as we know, there is no reliable method 
that can handle the very large configuration space arising 
in such problems. The HBO approach allows one to use 
a significantly reduced configuration space and is likely 
to be the only approach that is capable of solving such 
four-body problems for large mass ratios. 



VII. CONCLUDING REMARKS 

Concluding the paper, we would like to emphasize sev- 
eral issues important for future studies. First of all, mix- 
tures of light fermionic atoms with heavy fermions or 
bosons are promising for the observation of the Efimov 
effect. The search for Efimov trimer states was in the 
agenda of physicists for more than 30 years, mostly in 
nuclear and atomic physics [l^, [l^, [l3| • The main reason 
is a remarkable discrete scaling symmetry of the system: 
dimensionless three-body observables are invariant under 
the transformation a A" a, where n is any integer and 
the value of A is fixed by Eq. For example, invari- 

ant are the ratios of binding energies to eg or the ratio 
of the atom-dimer scattering length to a. However, the 



evidence for Efimov states was obtained only recently in 
Innsbruck experiments with cesium atoms through the 
measurements of the rate of three-body recombination 
as a function of a [3) [3 • An important aspect is that 
the factor A for a three-body system of one light and 
two heavy atoms can be relatively small, for example 
A w 5 for the i^^Cs-issCs-^Li system instead of 22.7 for 
the Cs3 trimer. This can allow the observation of several 
oscillations of the log-periodic dependence of three-body 
observables in the same window for the values of the two- 
body scattering length a. 

We should mention here that a direct observation of 
bound Efimov trimers can be rather difficult because of 
their short intrinsic lifetime. Indeed, considering a trimer 
with the energy |e,y| « 10/iK and assuming an optimisti- 
cally small ry* ~ 10~^, we obtain the trimer lifetime of 
the order of 100 /xs. Therefore, one should invent indi- 
rect ways of identifying the existence of Efimov trimers. 
The key idea that we bring in with this paper is that 
the observation of the Efimov effect can be made in a 
gas of dimers. From our results for the dimer-dimer in- 
elastic collisional rate (trimer formation) in Figs. [2] and 
[3] one can extract information on the binding energies 
of trimer states, their lifetimes, and also on the value 
of the three-body parameter. The dependence a{a) im- 
mediately gives the parameter ry* and, hence, the trimer 
lifetime r. The positions of the minima of the curve a{a) 
indicate the values of the two-body scattering length a 
at which one has a trimer state with the binding energy 
2eo. 

One can think of obtaining long-lived bound light- 
heavy-heavy trimer states in an optical lattice acting on 
heavy atoms and increasing their effective mass 27] . The 
heavy atoms in these states will be localized in different 
lattice sites and the light atom will be delocalized be- 
tween them. For example, one can consider the *°K-^°K- 
^Li system which does not exhibit the Efimov effect in 
free space and will show the presence of Efimov trimers in 
an optical lattice under an increase of the K mass by more 
than a factor of 2. The intrinsic lifetime of such trimers 
will be long as their decay into a ^Li-''°K (or ^"K-^"K) 
deeply bound molecule and a free atom requires two 
atoms to be on the same lattice site, which is strongly 
suppressed by the Pauli principle. 
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APPENDIX 

Explicit expressions for the operators L and L' in the 
case of s-wave scattering [/(R) = /o(^)] and p-wave scat- 



tering [/(R) = /i (i?) cos 6'r,] are obtained by performing 
the angular integration in Eqs. ([55)1 and ([55)1 . The cal- 
culation is straightforward and here we present only the 
resulting expressions. 



In the case of s-wave scattering we get: 



q sm a cos ( 



27r2 



P / dR 



,R 
B 



Ki {q\R-R'\/smt 



R-R 



{R ^ -R) 



(/o(i?)-/o(i?')) 



+ 



1 



cos 26* 



Ki {qVR'^ + R'^- ^RR' cos 261/ sin ( 
VR^ + R^ - 2RR' cos 29 



(R' ^ -R') 



MR') 



(50) 



Here q = -y/ AI {cq — E) and in each of the square brackets the second expression is the same as the first one with R 
replaced by — i?'. 

The action of the operator L' on a spherically symmetric function reads: 

07 I 1 POO 

L'fo{R)= Y.\-^r.AR)xUR')K^UR.R')-xlAR)xl:i^^^^ (51) 

1 = 1,3,,... n 

where Xn,i{R) is the radial part of the eigenfunction of Eq. (j20p with orbital angular momentum I and radial quantum 
number n. The quantities x° ;(i?) and k° ^ correspond to solutions of Eq. ([^D]) in which e+{R) is substituted by eo- 



The kernel K^^ is defined as 



K^\Ri,R2) - 2/,+i/2(«i?</2)i^i+i/2(^i?>/2)//Rli?^ 



(52) 



where the functions I^i{z) and K^z) are the modified Bessel functions of the first and second kind, and i?< {R>) is 
the smallest (largest) of the coordinates Ri,2- For = €n,i < E, i.e. for negative imaginary Kn,i given by Eq. ((M)) . 
we use the rules K{—iz) = (7ri'^+^/2)7?^(z) and Ifj,{—iz) — Jf^{z). 

In the p-wave case we define the radial operator by the relation L[/i(i?) cos^r] = cos0R,L'"^/i(i?). The 

explicit form is given by 



L'='hiR) 



q sm a cos f 



27r2 



P / dR 



,R 
R 



Ki{q\R- R'l/sint 
\R-R\ 



- {R ^ -R) 



MR) 



Koiq\R- R'\/sine) Ko{q\/R^ + R'^ - 2 RR cos 29/ sine) 



qRR 



cos2 29qRR 



+ (R ^ -R) 



MR') 



Ki{q\R- R\/sm9) Ki {qVR^ + R^ - 2RR' cos 26*/ sin ( 



(R ^ -R) 



MR') , (53) 



\R-R\ cos 29 A R^ + i?'2 - 2RR' cos 29 

and, accordingly, for the radial part of the operator L' we have 

{L'y=^MR)^^ E / T.^Xn,i{R)xUR')Kt,{R,R')~x%{R)xl:i{R')Kl^lS^^ (54) 



(=0,2,... 

where the kernel K^^^^ is given in terms of K^'''' from Eq 

RW ^ (i?, R) = {l + l)K(J+p{R, R) + iMJ-piR, R') 

I 



(55) 
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